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. . Abstract 

O \ We present new measurements of the thermal conductivity {k) of UPts down 
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to very low temperatures (16mK) and under magnetic fields (up to 4 T) which 
cover all the superconducting phases of UPts. The measurements in zero field 



> 

^^ are compared with recent theoretical predictions for the thermal conductivity. 
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^^ \ which is dominated by impurity states at the lowest temperatures studied. 

The measurements under magnetic field at low temperatures are surprising 

a _ 

C^ I since they don't show the expected low field square root dependence, k oc VB. 

'^ ' The discontinuity of dn/dT at Tc changes drastically when passing from the 
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high field low temperature C phase to the low field high temperature A phase 
: this is related to the change of the symmetry of the superconducting order 



H ' parameter when crossing the A— > C phase transition. 
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I. INTRODUCTION 

UPta is widely believed to be an unconventional superconductor. The existence of sev- 
eral different superconducting phases in the B- P-T phase diagranim is experimentally well 
confirmed, and can only be explained by a change of the symmetry of the order parameter 
(OP) when passing from one superconducting phase to another (for reviews, see e.g. Ref. 
|]I|]). This, in turn, can only be understood if additional symmetries, other than gauge sym- 
metry, are broken at the superconducting transition. The hexagonal crystal structure of 
this compound restricts the possible OP to 8 one-dimensional (so called A,B) and 4 two- 
dimensional (E) irreducible representations of the crystal point group, which are further 
classified according to the parity of the superconducting wave functions (index u for odd 
and g for even parity) [0]. In some of the irreducible representations the order parameter 
changes sign under some of the point group symmetries, implying points or lines of zeros of 
the gap on the Fermi surface. Much theoretical and experimental work has been done in 
order to try to find the symmetry of the OP in this compound, but up to now, no definite 
consensus emerges. The thermal conductivity is sensitive to the position and type of the 
nodes (line or point nodes) of the superconducting gap, which itself depends on the sym- 
metry of the OP. Previous work (See Refs. 0J^J^-0| for experimental and Refs. p|-p!9|] for 



theoretical work) has already demonstrated that the thermal conductivity (/t) in this com- 
pound is of electronic origin at very low temperature, and thus sensitive to the electronic 
excitations in the superconducting state. Compared to the specific heat C, k is a directional 
probe, and it seems |^,^ that the anomaly present in the specific heat below lOOmK does 
not influence k. Thus, a directional study of k down to very low temperatures is expected to 
give much information about the gap structure of UPta. Note that the first measurements 
0jl^ had already stressed the necessity of choosing a model for the OP yielding a hybrid 
gap with a line node in the basal plane and point node along the c axis. Measurements on 
far improved samples to lower temperatures are now expected to distinguish between such 
models. The discussion of our zero field data will be limited, as in all recent theoretical 



T5|-|19[ and experimental 0-§| work, to two plausible symmetries for the OP that appear 
consistent with the observed properties of superconductivity in UPta : the Eig and E2U 
representations of the order parameter which both lead to a hybrid gap, with a line node 
in the basal plane and nodes along the c axis where the gap vanishes respectively linearly 
(Eig) or quadratically (E2„) with the polar angle 9. Up to now, no experimental distinction 
could be made between the two models on the basis of the thermal conductivity. Since the 
work of K. Behnia et al. [^j measurements (Lussier et al. for T>100mK) emphasizing the 
appearance of an anisotropy of the thermal conductivity below Tc concluded that the gap 
was of the Eig type, and further measurements by the same group to lower temperatures 
(T>50mK) IP] were interpreted to favour E2U. We measured the thermal conductivity 
down to 35 mK, and showed that although some features of our data, like the power law 
behaviours, seemed to favour the E2U model, others did not, like the temperature depen- 



dence of the anisotropy. In fact it was shown [|T^-|T9| on the basis of the data of Lussier et 
al. PI that the existing experimental data could not differentiate between the two models. 
An important point is that for the lowest temperatures, which should be best suited to 
distinguish between the two gap structures, theory predicts that the pair breaking effect of 
impurities and defects produces a band of impurity bound states which restores a behaviour 
of the thermal conductivity comparable to that of the normal state. This regime was never 
observed up to now, as the measurements were not extended to low enough values of T/T^. 
Our data down to 16 niK show the onset of this type of behaviour, which is accompanied 
by drastic changes on the anisotropy of k below 30 mK. In addition, the experimental res- 
olution has been considerably improved, thus our data give a much better basis for a more 
detailed discussion of the thermal conductivity. We discuss extensively our zero field data, 
and also present new measurements under magnetic fields. Measurements under field in the 
mixed state give an additional means of probing the gap structure P,^pl|-[26[1 . Indeed, the 
possibility of a hybrid gap such as that in the Ei^ and E2U models was proposed by Yin and 
Maki [^ on the basis of the mixed state thermal conductivity data of Behnia et al. 



Our data are more systematic and show new features of k in the mixed state, arising both 



at very low temperatures and near the tetracritical point, where the three superconducting 
phases of UPta coalesce. 

II. EXPERIMENTAL 

The crystals used for the measurements were cut at adjacent sites of the same mother 
crystal and followed exactly the same heat treatment. Their critical temperature Tc is the 
same, their specific heat shows [^ clearly the double superconducting transition and their 
residual resistivities are among the lowest ever reported (po = 0.5AfiQcm for the current 
j//b, and po = 0.17fiQcm for j//c ||^). The samples were heated by a PtW alloy strain 
gauge heater, and the temperature was measured by two matched Matshushita 68^2 carbon 
thermometers. The accuracy of the measurements has been significantly improved. We 
could achieve a 1% resolution at 500 mK, 3% at 50 mK and 10% at 16mK for the measured 
anisotropy Hc/i^b- By comparison, the authors of Ref. [Q give error bars of 1% at 150 mK 
and 5% at 50 mK on n, whereas for the anisotropy Hc/i^b, the error bars are estimated to 
be 15% above 50 mK |18[. We improved our experimental method (described in Ref. 0) in 
the following points: (i) better electronics, and extensive time averaging to reduce noise. As 
a consequence, a continous run from 16mK to IK fully computer controlled took 3 days of 
measuring time, (ii) The long time averaging was possible due to the very high temperature 
stability of our dilution refrigerator, and its low temperature minimum : 5 mK. (iii) In 
this dilution refrigerator we could calibrate our thermometers during the same experiment 
against a Ge thermometer, a paramagnetic CMN salt and a ^He vapor pressure thermometer. 
The lowest measuring temperature for our samples was limited to 15 mK by their self heating 
due to the U radioactivity. The problem of the magnetoresistance of the thermometers (the 
improved sensitivity is payed for by important changes of the calibration under magnetic 
field) was solved by making a new calibration against a standard Ge thermometer in the zero 
field compensated region of the magnet at every field. The experimental error bars become 
more important when increasing the magnetic field, and will be shown in the figures when 



necessary. In high magnetic fields (3T), and for j//c, the thermal conductivity of UPts is 
extremely large : 3 orders of magnitude larger than at 16mK and zero field, due to the low 
value of the residual resistivity po for j//c (po =0.17 pficm). This, and the poor geometric 
factor of our sample for j//c (roughly 2 mm length , compared with 6 mm for the sample 
measured for j//b) conspired to make the measurements at 3 T and very low temperatures 
extremely difficult. Therefore, at 3 T and for j//c we only measured above lOOmK. Taking 
another sample with better geometry would have solved this problem, but our interest was 
focused on the superconducting phase of UPts. More than 60 k(T) runs under different 
magnetic fields were done in order to scan the B-T phase diagramm. All magnetic fields 
where applied along the same axis as the heat current (B//j//b and B//j//c). 

III. NORMAL PHASE 

The normal phase behaviour of the thermal conductivity of UPts shows some features 
important for the understanding of the superconducting phase. In Figure 1 we have plotted 
K,{T) for j//b and j//c, at T and at 3 T where the normal phase is recovered. The dominant 
scattering mechanisms of the electrons are the same as in the electrical conductivity : (i) 
elastic scattering of electrons by impurities and (ii) inelastic scattering between electrons. 
The first scattering mechanism (i) leads to a linear temperature dependence of the thermal 
conductivity and is related to the temperature independent part of the resistivity by the 
Wiedemann- Franz law : k/T = Lq/pq [Lq = 2.44 10~^W/K'^cm). As seen in Fig. 1 this 
mechanism dominates at very low temperatures T<100mK, below which n/T only varies 
by 10% as T— >0K. The ratio of the thermal conductivities extrapolated down to T=OK 
along the b and the c axis is 3. This ratio is usually said 0,|,^^-|T8| to be related to the 
anisotropy of the Fermi velocities, rather than to the anisotropy of the scattering rate r, 
which is assumed to be isotropic. Note that this point is not experimentally demonstrated. 
Nevertheless we will also take r to be isotropic, in order to compare our data with theory. 
The second mechanism (ii) leads to the well known Fermi liquid T^ behaviour of the electrical 



resisitivity and is of growing importance with increasing temperature. Deviations from the 
simple Wiedemann-Franz law are then expected : one expects again n/T = L/{pq + AT"^) 
but now with a temperature dependent Lorentz number L(T)< Lq. 

In order to verify this point, we also measured the electrical resistivity of our samples at 
3 T. In high purity samples, uJcT {uJc being the cyclotron frequency) becomes larger than 1 at 
low temperatures. This complicates the temperature dependence of the resistivity at 3 T, so 
that we did not get the expected simple Fermi liquid behaviour p = po + AT"^. Such effects 
have been already observed in UPta [^^. Nevertheless we can roughly estimate that the 
Lorentz number varies from 0.95 Lq at T=OK to 0.8 Lq at T=600 mK. We note that the 
positive magnetoresistance of UPta (po varies for j//b (j//c) by roughly 20% (10%) between 
and 3 T PJ2^) leads one to expect that k decreases with B, as was already observed by 
Behnia et al. but contrary to what was shown by Lussier et al. P], although these authors 
might have corrected qualitatively for the magnetoresistance, assuming the validity of the 
Wiedemann-Franz law in all the measured temperature range. 

IV. SUPERCONDUCTING PHASE 

The thermal conductivity in zero field is shown in Fig. 1. No significant change is observed 
at Tc. It is important to realize that up to now quantitative comparison with theory in 
the superconducting phase can only be done easily when elastic scattering with impurities 
prevails. As shown in Fig. 1, k at Tc is roughly half the value of the zero temperature 
extrapolation of n/T in the normal state. This means that at Tc both mechanisms (i) and 
(ii) are of equal importance. When entering the superconducting state, k rapidly drops as 
the superconducting gap opens over the Fermi surface. Therefore, the inelastic scattering 
rate in the superconducting state should diminish more rapidly than in the normal state. 
A qualitative treatment of the inelastic scattering rate has been proposed by Fledderjohann 
and Hirschfeld |jT^, extrapolating the normal phase scattering rate taken to be proportional 
to I/Ipq+AT"^). Norman and Hirschfeld |TB[ also made an extensive analysis of Fermi surface 
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effects on k, including inelastic scattering, and demonstrated the difficulty of fitting k in 
this temperature range. As suggested in Ref. [0, they tried to fit at the same time k(T) 
and the anisotropy Kc/z^fe, shown in Fig. 2. This quantity has little temperature dependence 
in the normal phase. As shown in the inset of Fig. 2, our enhanced resolution permits us to 
observe a kink in nd i^h against temperature at roughly 550 mK, which should correspond to 
the onset of superconductivity. Near Tc (for T<Tc), Kc/ i^b is independent of temperature, 
as shown in Ref. p, but contrary to what was observed by Lussier et al. [0,|1. Below 0.8T~ 
(T^=480 mK) k^ Hh starts to decrease and seems to tend towards a finite zero temperature 
extrapolation. However below 30 mK (0.07T~) a stronger dependence on temperature sets 
in. Norman and Hirschfeld calculated nd k\, using many different functions to describe 
the OP within the Eig and E2n irreducible representations, and one of the curves they 
calculate (Fig. 8 of Ref. |T^) fits our measurement of the anisotropy Kc/k;, well, however 
their calculated temperature dependence of k^ and Kc do not correspond with our data. 
Nevertheless, we note that a temperature independent anisotropy i^d ^h near T^ seems to 
be only reproduced within E2m in their calculations. At the lowest temperatures (T— i^OK) 
the thermal excitations which drive the heat current have k vectors corresponding only to 
the Fermi surface regions near the nodes of the gap, and inelastic scattering can be ignored, 
so that the interpretation of k is expected to be simpler. We now discuss impurity effects 
which turn out to be important in this limit. 

In non conventional superconductors with nodes of the gap, the pair breaking effect of 
even non magnetic impurities or defects leads to a finite density of states at the Fermi level, 
due to virtual bound states forming a band of width 7 (see Refs. [p|- [TT|jr3| - [18| , |3TH 3^). This 



band is easier to observe when the impurities are in the unitarity scattering limit (scattering 
phase shift 5 = vr/2), and in heavy fermion superconductors, impurity scattering is expected 
to be in this limit 0. Thus, in UPts, the thermodynamics is expected to be governed by 
the propeties of this band for T< 7. 

Before discussing the peculiar properties of this band which leads (see below) to a finite 
KJT for T— i>OK, we focus on the power law behaviour expected for intermediate tempera- 



tures: 'J < T <^ Tc. In this temperature region, the impurity induced band can be neglected, 
as well as inelastic scattering effects. In the case of resonant scattering of the thermal ex- 
citations by impurities, simple power laws are expected, controled by the gap nodes. Our 
measurements of n/T below 80 mK are shown in Fig. 3 as a function of T^. We can observe 
two regimes : (i) between 80 mK and 30 mK a T^ power law behaviour of k, for both axis, as 
expected for a hybrid gap with resonant scattering and (ii) below 30 niK the effects of the 
impurity induced band of quasiparticles leading to a finite value of k/T for T— i>OK. Note that 
the extrapolation of the data for T>30mK gives a negative zero temperature extrapolation 
for j//c (zero for j//b), which implies that the anisotropy Kc/i^b still depends on temperature 
down to 30 mK (see Fig. 2). In an earlier paper we argued that observing the same power 
laws for both axis is in favour of the E2u model, as the quadratic angular dependence of 
the nodes near the c axis leads to a high density of thermal excitations comparable to the 



density in the basal plane [0. Nevertheless, the theoretical analysis of the authors of Refs. 
p!5| , p!6|JT8| , p!9[] demonstrated that this observation could also be consistent with an Eig state. 
Barash and Svidzinsky []TU[ examined the expected power laws in detail, and they find to 
leading order in T: 

- for j//b (heat current in the basal plane ): Kh = T^ {ab + bbln^ (T / Aq) +Chln{T/ Aq)) for 
Eig and E2u; 

- for j//c (heat current perpendicular to the basal plane): Kc = T^(ac + hcln^iT / Aq) + 
cJniT / Aq)) for E2u and Kc = acT^ for Ei^; 

with coefficients a, b and c which depend on the parametrization of the gap near the 
nodes, and even on contributions of the order parameter from regions of the Fermi surface 



far from the nodes [0. The power laws expected on Kc for Eig and E2u differ mainly by 
a logarithmic factor (1 + —Iii?{T/Aq) + —ln{T/Ao)). As seen in Figure 3, both Kf, and Kc 
do not follow a pure power law, but need a negative zero temperature extrapolation to be 
fitted. The formulas given above (for E2u) fit our data well, but with negative coefficients 
bf, and be, which may not have a clear physical meaning. Therefore, on the basis of power 
law behaviours we are not able to discriminate unambiguously between Ei^ and E2„. 
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The lowest temperature regime (T< 7) might still be used for the identification of the 
order parameter. Indeed theoretical work has shown |T5|JT6| , p!8[| that the existence of the 
impurity induced band of quasiparticles does not wash out completely the sensitivity of the 
thermal conductivity to the choice of the OP. We will follow closely the discussion of Graf 
et al. |T^. They predict for T < 7 : Hi/T = a^ + /?jT^ with ctj and /3i which depend on the 
choice of the OP and of the direction of the heat current : 

- For i=b (heat current in the basal plane) ki, is the same for Eig and E2U, as both 
predict a line node in the basal plane. Moreover, a^ is universal, i.e. independent of the 



impurity concentration rii although it depends on the form of the gap near the node, and 
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- For i=c, ttc oc 7 and j3c = 2.5^ac for Eig, but for E2„, ac oc a^ is universal and 

The inset of Fig. 3 shows k/T as a function of T^ at the lowest temperatures (T<35 
mK) together with the fits to the predicted formulas. For the thermal conductivity in 
the basal plane, the fit gives at, = O.lSmW/K'^cm, of the same order of magnitude as 
ab ~ lmW/{K'^cm) estimated in Ref. [|K|. It also yields (3b/c(b = 4.2 lO^/K"^, therefore 
7 ^ llniK. Note that with a less pure sample we should have been able to observe the 
expected law over a larger range of temperatures. Assuming that the impurity concentration 
Uj is the same for both measured crystals, and having fixed a^ and jSb/ab only Oc can be varied 
as a free parameter in order to reproduce the measured curve for j//c [|15[. For E2U the best 
agreement is obtained with a^ = 0.2mW/K'^cm and for Ei^ with a^ = O.llmW/K'^cm. As 
shown in Fig. 3, the fit to Eig is better, but only in a small temperature range {16mK < T < 
25mK). Nevertheless, the assumption that 7 is the same for both measured samples may 
not be justified. We recall that the assumption of an isotropic scattering rate may not hold, 
and that even if both samples come from the same mother crystal and had exactly the same 
heat treatment, differences in the impurity (or defect) concentration can never he excluded. 
Indeed, although the onset of a finite resistivity occurs at the same temperature for both 
samples, the width of the superconducting transition measured by resistivity is 17mK for 
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j//b, but only 8mK for j//c (10%-90% criterion). The broadening of the superconducting 
transition could be correlated to the sensitivity to stress produced in the cutting of the 
material despite the subsequent annealing. Note that our measurements under magnetic 
fields show a clear anomaly in k, at Tc. As shown in the inset of Fig. 4 this anomaly 
coincides with the onset of a finite resistivity, thus with the lowest part of the resistive 
transition. Therefore, differences in the transition width, although they show that the 
samples are different, do not give much information about the defect concentration in the 
bulk of the material. We note that in order to explain our data within E2u, the impurity 
concentration Uj oc 7^ ||15[ should be 2.5 times lower for j//c than for j//b. One way of 
solving this problem is to measure less pure samples (in order to observe the expected laws 
in a larger range of temperatures), again for both axis and down to the lowest possible 
temperatures. This will also give an experimental test of the prediction of a universal value 
of K,/T{T -^ OK) for the heat current in the basal plane. Note that the upper critical field 
anisotropy (6^2 is slightly Pauli limited for B//c but not for B±c) is an important argument 
in favour of an odd parity OP in UPta (see Ref. |^^). It will be very difficult to explain Bc2 
within Eig. 

One important point not understood at present is the discrepancy between the value of 7 
estimated from the normal phase data (with the impurity scattering rate Tq/Tc ~ 5 10^^ we 
obtain 7 ^ 50mK] see Ref. [|^ for the formulas), which is 3 times larger than the one found 



from our low temperature fit. Indeed, as pointed out by Norman and Hirschfeld |T8[ the data 
at higher temperatures (7 ^ T < Tc) seem to be difficult to fit with such a low impurity 
scattering rate. This may suggest that the description of the impurities and defects in UPta 
only in terms of s-wave scattering in the unitarity limit is too simplified for the discussion 
of subtle very low temperature effects. It also stresses the lack of experimental knowledge 
on the nature of the dominant scattering centers in this compound. 
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V. MIXED PHASE 

Figure 4 shows an example for our temperature scans at different magnetic fields. The 
first striking feature is that under magnetic fields an anomaly appears, which was not visible 
in zero field. More precisely, there is a jump in the derivative of /t : d/t/dT at Tc. This 
is predicted to happen already in conventional s-wave superconductors, and the interesting 
point is that, from the superconducting side, d/t/dT is related to the order parameter. At 
intermediate temperatures (100 mK<T<Tc), k, first decreases at low fields as a function 
of B, and then increases up to its value in the normal phase (this is seen more clearly in 
the Fig. 8, where we plot /t as a function of B). At very low temperatures (T<50mK), 
no decrease at low fields is observed and k, increases continuously with B. Qualitatively 
speaking, there are two competing effects : an enhancement of the thermal conductivity 
with the magnetic field towards the normal phase value, and the scattering of the thermal 
excitations present outside the vortex cores by the vortices which may tend to diminish k 
(compared to the zero field value). The understanding of the thermal conductivity of clean 
superconductors in the mixed state is difficult due to the complex interplay between these 
effects, and new theoretical calculations are needed in order to understand ^(B). Nevertheless 
some conclusions can be drawn. 

All of these effects strongly depend on the relative direction of the magnetic field and 
the heat current. Our choice was to measure k with the field always parallel to the heat 
current, along the c axis and in the basal plane, in order to make the interpretation easier. 
Also, the magnetic field was always changed in the normal phase (field cooled). We checked 
that only a small irreversibility exists at very low fields (B<0.03T). 

Before discussing the data, we note that we recover the phase diagram of UPts when 
plotting Tc obtained by the anomaly of /t as a function of the magnetic field (see Fig. 5). The 
B^C line could also be observed due to anomalies shown in Fig. 8 on the field dependence 
of K, but no clear signature could be found of the A^B phase transition. We will first 
discuss the thermal conductivity near Tc(B), then focus on k{B) as T^O K, and finally 
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discuss k{B) at intermediate temperatures. 

In the literature (see Refs. [P,^ p5| , p6| for the discussion of UPta), it is generally the jump 



in dn/dB which is discussed. We have measured dn/dT because our resolution for temper- 
ature scans was much better than for field scans, but the main conclusions should remain 
the same. In Fig. 6 we show the measured dn/dT at T,, as a function of the magnetic field, 
normalized to the normal state values : {dug/ dT — dKn/ dT) / {K,n/Tc{B)) , where the subscript 
s means approaching Tc from the superconducting state, and n from the normal state. As 
already mentioned pSJ^ and shown by the measurements of Behnia et al. |P,^ in UPts, the 



slope dus/dB, and therefore also dK^/dT, depends on the direction of current and field, and 
is related to the gap structure. Therefore, if the topology of the gap changes when passing 
from the high temperature low field A phase to the low temperature high field C phase, a 
jump of dn/dT should be observed at the A— >C phase transition. Indeed, as shown in Fig. 6, 
we clearly observe this jump for both orientations of the field at the A^C phase transition 
(0.5 T for B//b and 1 T for B//c). If no phase transition took place, the overall behaviour 
of the measured quantity would be parabolic like, with {dn^/dT — dun/ dT) / {i^n/Tc{B)) = 
at B =0 and aX B = Bf^{T = 0), with a maximum at an intermediate field. We expect this 
parabola to depend on the topology of the superconducting gap. So it might be possible to 
draw more conclusions about the gap structure of the A and the C phases, from which little 
is known, by a careful theoretical analysis of these data. 

At very low temperatures (T<50mK) we observe a continuous increase (Fig. 7) of n{B) 
up to its normal phase behaviour, which we will discuss now. First, we recall the situation 
in clean s-wave superconductors. Within a naive model, one might think that in the limit 
T— i>OK the contribution to the thermal conductivity coming from the electronic excitations 
within the vortex cores would dominate n, as this type of excitation dominates the density of 
states (in s-wave superconductors). One would then expect that k scales roughly as -^^n, 
but this turns out not to be correct. Indeed, the low group velocity of the excitations within 
the vortex cores leads to a small contribution of this type of excitations to the thermal 
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conductivity. Vinen et al. [ESl estimate Kvortex cores ~ 10 ^-^/tn for pure Nb. We observe in 
UPta an increase of k with the magnetic field which scales roughly as -^i^n for B//j//b and 



as ^^i^n for B//j//c at low fields (Fig. 7). This is a surprising result which is not explained 
with the existing models. It might be related to the peculiar magnetic field behaviour of 
the density of states of unconventional superconductors with a line of nodes of the gap. 
Indeed, it has been predicted that in such a case, the density of states would be dominated 
by the contribution of the thermal excitations of the superconducting phase outside the 
vortex cores, leading to N(Ei7) oc v5 instead of N(Eir) oc B when the main excitations are 



those in the vortex cores. Note that this is valid at low fields for Bd < B <^ Bc2 ||30| , p^ 
Again, in a naive model, one would expect that ^(B) as T-^OK follows the density of states, 
but no signature of the expected square root behaviour due to the line of nodes predicted 
within Eig and E2u is found in our data. Nevertheless, it is already known that in the case 
of clean s-wave superconductors, the thermal conductivity does not in general follow the 
density of states. We note also that the thermal excitations outside the vortex cores are 
scattered by the vortices (see below), and this may considerably infiuence |;(B) also in the 
limit T— >0K. A rough square root behaviour was found in the specific heat C by Ramirez 
et al. 1^, although at suspiciously high temperatures (150 mK; Tc/3), and large fields (up 
to Bc2/3). These measurements are therefore infiuenced by the field dependence of T^ and 
by a huge low temperature anomaly present in C [§3'§M1 • "^^^ origin of this anomaly, which 
appears already at 100 mK in zero field is not understood at present, but it is for sure not 



connected with the superconducting state [^. So the measurements of Ramirez et al. ||2^ 



are clearly not simply related with the field dependence of the density of states at T^OK. 
The advantage of the thermal conductivity is that it seems not to be infiuenced by this 
anomaly. Its origin may therefore be related to localized modes which do not carry heat . 
In summary, it is clear that k{B) as T^OK has a peculiar behaviour in UPts not explained 
by the existing models. One must explain not only the linear behaviour of k, with B at low 
fields, but also the anisotropy found for the constants of proportionality {k, ^ B/Bc2 for 
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B//j//b and K ^ \B/B,2 for B//j//c). 

As /t(-B) scales roughly as ^^f^n for B//j//c, ^(5) cannot remain linear in B up to Bc2. 
Indeed, at roughly 0.8Bc2, i^c strongly deviates from the linear behaviour in order to reach 
the normal state value at Bc2, whereas k^ is roughly linear in the whole field range up to Bc2. 
This behaviour is observed for the extrapolation T— >0K, as well as at 100 mK, as shown 
in Fig. 8. Note that Behnia et al. P] already observed a large value for d/t/dB with B//c 
when the upper critical field is Pauli limited along the c axis |^. This may influence k(B) 
near Bc2- This measurement, and the ones of Behnia et al. (Refs. P,^ in UPts and Ref. ||2^ 
in URu2Si2, another heavy fermion superconductor with a Pauli limited upper critical field 
along one crystallographic axis) are to our knowledge the only existing thermal conductivity 
measurements in the mixed phase of superconductors with a Pauli limited upper critical 
field. We could not find theoretical predictions in this case for k,{B) near Bc2- 

We now focus on the magnetic field behaviour of k between 100 niK and Tc, shown 
in Fig. 8 and obtained by plotting our temperature scans as a function of B. As pointed 
out above, the scattering of thermal excitations by vortices leads to a diminution of k{B) 
visible at low fields. This type of scattering was indeed reported experimentally in s- wave 
superconductors (see e.g. Refs. PU| , P7| ), and in UPts by Behnia et al. [^ and Huxley et 



al. 1^. In a simple model for low fields where vortices can be treated as independent 
scattering centers, the additional thermal resistivity W=l/fi: (Fig. 9) coming from scattering 
of the superconducting thermal excitations by the vortices is expected to be proportional 
to the magnetic field (see e.g. Refs. ||20| , |27|j36| ) . In s-wave superconductors, the law W = 
W(0)(l+^i?) (/o is the mean free path in the superconducting phase, $o the flux quantum, 
and a the effective scattering width of the vortices) is observed at low fields and in a large 
range of temperatures. Our data do not show such a behaviour. To understand this, the field 
dependence of n for T— i^OK has first to be explained. The fact that we do not observe the 
classical behaviour may indeed be related to the unconventional nature of superconductivity 
in UPts. 
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VI. CONCLUSION 

We have shown that the thermal conductivity shows a rich behaviour in the different 
superconducting phases of UPta. Experimentally, we could reach low enough temperatures 
to enter the regime where the thermal conductivity is governed by impurity bound states. 
It is unclear if it is now possible to differentiate between Ei^ or E2„ on the basis of this data, 
we can however parametrize the gap nodes for both models. It appears that measurements 
on samples of different purity down to the same lowest temperatures would be very useful. 
First, they could test the validity of the theoretical models, predicting for Eig as well as for 
E2U a universal value of k/T for T^OK in the basal plane. Second, they might distinguish 
between Ei^ and E2U '■ theory predicts for the heat current along the c axis a universal value 
of k/T for T^OK in the case of E2„, but a sample dependent value for Eig. Although our 
very low temperature data are slightly better explained within Ei^, we remind the reader 
that Eig has difficulties to explain other measurements on UPts such as the upper critical 



field [^ or the Knight shift [^ which both point towards an odd parity OP. 

We do not observe any coupling between the low temperature specific heat anomaly and 
the thermal conductivity at very low temperatures. This points to a localized nature of 
the modes that lead to this specific heat anomaly. As indicated by Fomin and Flouquet 
PP| it could be possible that the real magnetic ordering in UPts arises only at much lower 
temperatures (T ^ 20 mK) than indicated by the appearance of a neutron diffraction signal 
(at T ^ 5K). This signal could indeed be related to the onset of short range order as 
a consequence of large magnetic fluctuations, rather than to a real long range magnetic 
order |^. But it can not be excluded that a weak interplay between these local moments 
and the heavy electrons may influence the temperature dependence of k and complicate a 
quantitative fit. 

We have also presented measurements of the overall temperature and field behaviour of 
the thermal conductivity in the mixed state. The low field low temperature result shows 
the absence of the expected square root behaviour of the thermal conductivity, for both 
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measured axis. We have also seen that the particular field dependence of k,{B) near Bc2 
is very different along the b or c axis axis {dK(./dB ^ dnii/dB), which may be related to 
the Pauli limitation of the upper critical field along the c axis. Clearly, both points invite 
more theoretical calculations and might be deeply connected to the OP of UPts. In the 
intermediate field regime, both the vortex core contribution to k, and quasiparticle-vortex 
scattering are important, and the data are difficult to interpret without further theoretical 
work. Near Tc, the derivative of the thermal conductivity, directly related to the order 
parameter, shows a clear jump when crossing the A— ^C phase boundary. This is a clear 
signature of the change of symmetry of the OP at this phase transition. 
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FIGURES 
FIG. 1. The thermal conductivity of UPts below 1 K for the heat current along the b and along 

the c axis. The zero field data are shown together with the normal phase data, measured at 3 T. 

FIG. 2. The anisotropy Kc/k^ normalized to 1 at T^ = 480mK is shown as a function of the 
reduced temperature T/T~. The inset shows a kink of Kc/kj, at T^=550 mK. Note the low error 
bars compared with earlier work (Refs. PJ5l^,0]). 



FIG. 3. The very low temperature thermal conductivity of UPts (between 70mK and 16mK) 
as a function of T^, with the fits based on the predictions of Graf et al. |T^. The inset shows 
k/T between 34mK and 16mK. Below 30 mK, the properties of the impurity induced band of 
quasiparticles are observed (see text). 

FIG. 4. Some of our temperature scans at constant magnetic fields are shown in this figure. 
The inset shows the observed anomaly at Tc, compared to the resistive transition at 1.56 T. 

FIG. 5. Phase diagramm of UPta deduced from the thermal conductivity measurements. The 
A^N and C^N transitions are traced by the anomaly on the temperature scans which appear at 
Tc under magnetic field. The B^C phase transition is seen on the field dependence of k (see Fig. 
8). No clear sign of the A^B transition was seen in our measurements. Lines are guides to the 
eye. 

FIG. 6. The slopes obtained with our temperature scans under magnetic fields. Lines are guides 
to the eye : continuous lines in the A phase, and dashed dotted lines in the C phase. The error 
bars, which become more important at high fields due to the broadening of the transition are also 
shown. 

FIG. 7. The figure shows k/T{B) for T^OK at low fields. Note that surprisingly, Kb and Kc 
show a linear behaviour (see text). At very low fields, close to Bd a curvature of k{B) is observed, 
but in this region a non negligible irreversibility exist. 



20 



FIG. 8. The magnetic field dependence of k/T{B) at some representative temperatures together 
with the zero temperature extrapolation of n/T{B) (open circles). Note that, while Hb/T is roughly 
linear with the magnetic field, Kc/T shows a strong upturn near Bc2 up to the normal phase value. 

FIG. 9. The magnetic field dependence of the thermal resistivity W{B) = 1/k{B) at low fields 
normalized to the value at zero field Wq. Lines are guides to the eye. Note that our measurements, 
more extensive than in Refs. [§,||], show that W is not linear in B at low fields (see text). 



21 



This figure "Fig_l.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_2.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_3.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_4.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_5.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_6.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_7.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_8.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



This figure "Fig_9.GIF" is available in "GIF" format from: 



http://arXiv.org/ps/cond-mat/9702089vl 



